International Journal of Modern Physics A, Vol. 19, No. 5 (2004) 761-774 
© World Scientific Publishing Company 



INCREASE OF THE 

ELECTRON-POSITRON HAWKING RADIATION 
FROM SCHWARZSCHILD BLACK HOLES 
BY DIRAC MONOPOLES 



YU. P. GONCHAROV 
Theoretical Group, Experimental Physics Department, State Polytechnical University 
Sankt- Petersburg 195251, Russia 

N. E. FIRSOVA 
Institute for Problems of Mechanical Engineering 
Russian Academy of Sciences, Sankt- Petersburg 199178, Russia 



Received 18 November 2003 

An algorithm for numerical computation of the barrier transparency for the potentials 
surrounding Schwarzschild black holes is described for massive spinor particles. It is 
then applied to calculate the total (including electronic neutrino and the contributions of 
twisted field configurations connected with Dirac monopoles) luminosity for the electron- 
positron Hawking radiation from a Schwarzschild black hole with mass M = 10^5 g. It 
is found that the contribution due to monopoles can be of order 12 % of the total 
electron-positron luminosity. 
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1. Introduction 

The present paper is a natural continuation of our works ^'^ so, referring for 
more details on motivation of the given trend to the mentioned references, we shall 
here restrict ourselves to a short discussion. 

Although it has elapsed over 25 years since Hawking discovered the possibility 
for black holes to radiate quantum particles, the situation with quantitative study of 
the Hawking radiation up to recently should be considered as extremely unsatisfac- 
tory. If qualitative understanding of the process developed more or less successfully 
then attempts of quantitative calculations might be enumerated by fingers. Really 
one may speak only about works of Refs."^ and those of Refs.'* based on Page's 
ones. However, when analysing the mentioned papers one can note that central 
problem of numerical calculations for Hawking radiation - calculating the barrier 
transparency for the potentials surrounding black holes - was actually not resolved 
so it is impossible to extract some exact algorithm for numerical computation of 
the necessary quantities from the mentioned references. 

Another aspect of the problem in question is connected with that any isolated 
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black hole might possess the internal magnetic fields of the Dirac monopole types. 
The latter configurations should be connected with nontrivial topological proper- 
ties of black holes and could have an essential influence on quantum processes near 
black holes, for instance, on Hawking radiation. A number of examples of such 

configurations may be found in Rcfs.^ and references therein. Physically, the exis- 
tence of those configurations should be obliged to the natural presence of magnetic 
U(N)-monopoles (with A'' > 1) on black holes though the total (internal) mag- 
netic charge (abelian or nonabelian) of black hole remains equal to zero. One can 
consider that monopolcs reside in black holes as quantum objects withoiit having 
influence on the black hole metrics. They could reside in the form of monopole 
gas in which the process of permanent creation and annihilation of the virtual 
monopole-antimonopole pairs occurs so that the summed internal magnetic charge 
(i. e., related with topological properties) is equal to zero while the external one 
(not connected with topological properties) may difiier from zero. While existing 
the virtual monopole-antimonopole pair can interact with a particle and, by this, 
increasing the Hawking radiation (sec Refs.^ and references therein). 

In other words one may say that the black holes due to their nontrivial topolog- 
ical properties can actually carry the whole spectrum of topologically ineriuivalent 
configurations (TICs) for miscellaneous fields, in the first turn, complex scalar and 
spinor ones. The mentioned TICs can markedly modify the Hawking radiation 
from black holes. Physically, the existence of TICs should be obliged to the natural 
presence of magnetic U(N)-monopoles (with A > 1) on black holes and additional 
contributions to the Hawking radiation exist due to the additional scalar or spinor 
particles leaving the black hole because of the interaction with monopoles so the 
conforming radiation can be called the monopole Hawking radiation.^ 

Up to now, however, only influence of the TICs of complex scalar fleld on Hawk- 
ing radiation has been studied more or less (see Refs."*^ and references therein for 
more details). The description of TICs for spinors was obtained in Refs.^ but the 
detailed analysis of the TICs contribution to Hawking radiation requires knowledge 
of the conforming ^-matrices which regulate the spinor particle passing through the 
potential barrier surrounding black hole. Those S-matrices for fhc Schwarzschild 
(SW) black holes have only recently been explored in Ref.^ which allows us to ob- 
tain an algorithm to calculate the 5-matrix elements numerically since for physical 
results to be obtained one needs to apply the numerical methods. In the massless 
case the S'-matrices discussed are simpler to treat and the paper of Ref.^ contained 
a description of the algorithm needed and applied it to calculate the all configu- 
rations luminosity for massless spinor particles for a SW black hole. The results 
obtained can serve as an estimate, in the first turn, of the electron-positron Hawk- 
ing radiation and also for the neutrino one from the SW black holes. The more 
exact computation should take into account the particle masses and require more 
complicated algorithms for calculating the corresponding 5-matrices. 

The present paper contains a description of one such an algorithm and applies 
the latter to calculate the total luminosity for electrons, positrons, electronic neu- 
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trino and antineutrino for a SW black hole with mass M = lO'^^ g. Section 2 
is devoted to the common statement of the problem while Section 3 gives more 
detailed information about the relations needed for further calculations. Section 
4 estimates the correctness of computations and Section 5 presents the numerical 
results obtained. Finally, Section 6 is devoted to concluding remarks. 

We write down the black hole metric under discussion (using the ordinary set 
of local coordinates t, r, i?, (f) in the form 

ds^ = adt^ - a-^dr^ - r^{d-d^ + sin^ Mip"^) (1) 

with a = 1 — 2M/r and M is the black hole mass. 

Throughout the paper we employ the system of units with h = c = G = 1, 
unless explicitly stated otherwise. We shall denote L2{F) the set of the modulo 
square integrablc complex functions on any manifold F furnished with an integration 
measure while L^i^F) will be the n-fold direct product of L^i^F) endowed with the 
obvious scalar product. 

Finally, when computing we use the following value of the electron-positron mass 
(see Ref.^) 

Ato(e±) = 0.51099906 Mev . 



2. Preliminaries 

As was disscussed in Refs.^, TICs of a spinor field on black holes are conditioned 
by the availability of a countable number of the twisted spinor bundles over the 

X -topology underlying the 4D black hole physics. From a physical point 
of view the appearance of spinor twisted configurations is linked with the natural 
presence of Dirac monopoles that play the role of connections in the complex line 
bundles corresponding to the twisted spinor bimdles. Under the circumstances each 
TIC corresponds to sections of the corresponding spinor bundle E, which can be 
characterized by its Chern number n e Z (the set of integers). Using the fact that all 
the mentioned bundles can be trivilized over the chart of local coordinates (t, r, t?, </?) 
covering almost the whole manifold x §^ one can obtain a suitable Dirac equation 
on the given chart for TIC ^ with mass /Uq and Chern number n e Z that looks as 
follows 

Vr,^ = /io*, (2) 

with the twisted Dirac operator 2?„ = ij'^'V'^ and we can call (standard) spinors 
corresponding to n = untwisted while the rest of the spinors with n ^ should be 
referred to as twisted. Referring for details and for explicit form of Vn to Refs.^, it 
should be noted here that in i^lR^ x ) there is a basis from the solutions of (2) 
in the form 

^"~x/2^ VF2(r,a;,A)ai$Amv' ' 

where cri is the Pauli matrix, the 2D spinor = ^\m{^,<fi) = {^ixm,^2\m) is 
the eigenspinor of the twisted euclidean Dirac operator with Chern number n on 
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the unit sphere with the eigenvalue A = ±\/ {1 + 1)^ — while —I < m < I + 1, 
I > \n\. As for the functions Fi^2, they obey the system of equations 



V^drF, + + ^)Fi= iiflo - C)F2, 



(4) 



with c = uj/y/a and a of (1). The explicit form of the 2D spinor is inessential 
in the given paper and can be found in Refs.^. One can only notice here that they 
can be subject to the normalization condition at n fixed 



7r 27r 



b b 



and these spinors form an orthonormal basis in igl^^) ^"^y n S Z. 

By passing on to the Regge-Whcclcr variable = r + 2Mln(r/2M — 1) and 
by going to the quantities x = r^jM.y = r/M,k = wM, /x = /xqM, we shall have 
X = y + 21n(0.5t/ — 1), so that y{x) is given implicitly by the latter relation (i.e., 
— oo < x < oo, 2 < y < oo) with 

y' = dy/dx = l-2/y={y-2)/y = ao{x) (5) 

and the system (4) can be rewriten as follows 

{E[+ aiEi = biE2 , , . 

withiJi = Ei{x,k,X) = F+{Mx), F+{r*) = Fi[r{r*)],E2 = E2{x,k,X) = iF_{Mx),F. 
F2[r{r*)] and 

ai,2 = ± -\/ao , (7) 

y 

^1,2 = M\/ao T • (8) 



To evaluate luminosity of the Hawking radiation for spinor particles it is necessary 
to know the asymptotics of the functions £^1,2 at a; — > +00. As was shown in Ref.^, 
the latter asymptotics look as follows 

El ~iyAr^sii(fc,A)e*'=^^e''''"(2'=^^^ x ^ +00 , (9) 
E2 - ^fcTMSii(fc, A)e*'=''^e*^'"(2fe+r.)^ x ^ +00 (10) 

with 

where sii(k, A) is an clement of the S'-matrix connected with some scattering prob- 
lem for a Schrodinger-like equation (see Section 3). After this one can obtain (in 
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usual units) the luminosity L{n) with respect to the Hawking radiation for TIC 
with the Chern number n in the form (for more details see Refs.^'^) 

oo 

L(n) = A^^ 2(^ + 1)/ 

±A l=\n\ ^ 

oo 

r I I ^ ^ ' 

where A = (#)^^^ « 0.251455 • lO'^^ erg • s'^ ■ M^^ and M in g while /x = 

3.762426569^10-i^AtoM if jiQ in MeV and M in g. 

Luminosity L[n) can be interpreted, as usual, ^'"^ as an additional contribution to 
the Hawking radiation due to the additional spinor particles leaving the black hole 
because of the interaction with monopoles and, as already mentioned in Section 1, 
the conforming radiation can be called the monopole Hawking radiation.^ 

Under this situation, for the all configurations luminosity L of black hole with 
respect to the Hawking radiation concerning the spinor field to be obtained, one 
should sum up over all n, i. e. 

oo 

i = ^L(n) = L(0) + 2^L(n) (12) 

neZ n=l 

since L{—n) = L{n). 

It should be emphasized that in (11), generally speaking, sii(fc. A) ^ sii(/c, —A). 
Obviously, for neutrino there exists only L{0) since neutrino does not interact with 
monopoles (when neglecting a possible insignificantly small magnetic moment of 
neutrino). 

3. The scattering problem and description of algorithm 



3.1. The scattering problem 

It is evident that for numerical computation of all the above luminosities one 
needs to have some algorithm for calculating sii(fc. A). As was mentioned in Section 
2, sii(fc. A) is an element of the S-matrix connected with some scattering problem 
for a Schrodinger-like equation. The latter equation, as was shown in Ref.^, looks 
as follows 

u" + {k^ - ii^)u = qou (13) 

with potential 

A^ /oq 

qo{x, k, A, /x) = ^ % (x, k, A, /x) + {^/a^ - l)/x^ (14) 

y 
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while 

qix, k, A, = (l - ^) + (a + ^gi) - ^ (Iq, + ga) 
with 

_ mVoq _ IJ'Vaoi^/y - 2) 

As an example, in Figs. 1, 2 the numerically computed potentials qQ{x, k, A, fi) with 
A > and A < are shown. It should be emphasized that potential go of (14) 
satisfies the condition^ 

+ C30 

J \qo\dx < oo 

— oo 

only at /i = 0, so when /i 7^ the correct statement of the scattering problem for 
go should be modified. Besides one can notice that when A > 0, go (a;, A) > at any 
X e [—00, 00] while at A < 0, go(a;, A) can change the sign. 



Figure 1: Typical potential barrier for A > 0. 



Figure 2: Typical potential barrier for A < 0. 

In its turn, as was shown in Ref.^, the correct statement of the mentioned 
scattering problem for Eq. (13) consists in searching for two solutions u~^{x,k, X), 
u~{x, k, A) of Eq. (13) obeying the following conditions 

,,+ (^ . ^^_ / e^'" + Mk, A)e-'=- + 0(1), X ^ -00, 

-( , _ /■S22(A;,A)e-'''^ +0(1), x ^ -00, 

u ia;,fi;,A; - |^_^^^^^2ifc+a;) + S2i(A;,A)w;i^_i(-2ifc+a;) + o(l), x ^ +00, ^^^> 

where the functions 1 (±2) are related to the Whittaker functions (±2;) 
(concerning the latter ones see e. g. Ref.^) by the relation 



w±i0,d±z) = W^^0l{±z)e- 



K|3/2 



so that one can easily gain asymptotics (using the corresponding ones for Whittaker 
functions^) 



w. 



if^^i{-2ik+x) = e''^+-e*'^'"|2'=^^l[l + 0{\k+x\-% x ^ +00 , 
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w_i0^i{2ik+x) = e-'''^''e-'^^''\'^''^''\[l + 0{\k+x\-'^)], x ^ +00 . (16) 

We can see that there arises some S'-matrix with elements Sij,i, j = 1,2. As 
has been seen above, for calculating the Hawking radiation we need the coefficient 
sii(A:, A), consequently, we need to have an algorithm for numerical computation of 
it inasmuch as the latter cannot be evaluated in exact form. The given algorithm 
can be extracted from the results of Ref.^. 

3.2. Description of algorithm 

To be more precise 

sii{k, A) = 2ik/[f-{x, k, A), f+{x, k, A)] , (17) 

where [,] signifies the Wronskian of hmctions /~ , /+ , the so-called Jost type so- 
lutions of Eq. (13). In their turn, these functions and their derivatives obey the 
certain integral equations. Since the Wronskian does not depend on x one can take 
the following form of the mentioned integral equations 

xo 

f-{xo,k,X)=e-"'''° + ^ J sm[k{xo -t)]q-{t,k,X)f-{t,k,X)dt, (18) 

— CSO 

Xq 

(/-);(xo,/c,A) = -ifce-*'="''+ J cos[k{xo-t)]q-{t,k,X)f-it,k,X)dt, (19) 



and 



/+(a;o,A;, A) = u;^^ 1 (-2iA;+a;o)+ 
^ j lin[Wi^i {-2ik+xo)w_i0^i{2ik+t)]q+{t,k, X)f+{t,k,X)dt , 



k+ 

Xo 



(20) 



{f^yxixo,k,X) = -^Wi0^i{-2ik+xo)+ 



xo 

where 



dx 

+00 



^«^i/3,i (-2«fc+xo)w_i^_ 1 {2ik+t) 



q+{t,k,X)f+{t,k,X)dt, (21) 



q-{x, k, A) = ^^q + ao/x' , q+{x, k, A) = ^^q + 2^^^ (22) 

y y xy 

with q = q{x, k, A, /x) from (14). 

The potential q~ exponentially tends to zero when a; ^ — 00 and potential 

(7+ behaves as 0(x~^) as a; -l-oo. So one can notice that these potentials are 
integrable when x ^ — ooorx^ 00 respectively. The point x = xq should be 
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chosen from the considerations of the computational convenience. The relations 
(18)-(22) can be employed for numerical calculation of Su. It should be noted, 
however, that while Eqs. (18)-(19) are suitable for direct calculation this is not the 
case for Eqs. (20)-(21) — it is difficult to evaluate the functions iu{z) in (20)-(21) 
since we have yet no effective fast method to compute the Whittakcr functions W(z) 
at z required. In the present paper we, therefore, stick to the following strategy. 
We employ the asymptotic expressions of w-functions at a; — > +00 from Eq. (16) 
and replace Eqs. (20)-(21) by 



f+{xo,k,X) = e*(fc+^+''W2fe+^o)) + 



+00 



^ / Im[e'(*=^(^°-*)+'''"(^<'/*»]g+(t,fc,A)/+(t,fc,A)di, 



(23) 



+00 

/ 



Im 



(/+);(X0,fc, A) = ie'(*=+^o+/.ln(2fe+xo))(fc+ ^ + 



-Xq 



i{k+ (xo-t)+iJ,ln(xo/t)) 



q+{t,k,X)f+it,k,X)dt. (24) 



The latter equations are appropriate for numerical evaluation but under this ap- 
proximation, as numerical experiment shows, reliable calculations are only possible 
for fj, < 0.01. This upper limit corresponds to /Uo(e^) and M ~ 10^^ g and it is 
obtained when requiring that calculations lead to the smooth monotonous graphs 
for r(A:,±A) = |sii(fc,±A)p without any oscillations, breaks or fractures (see below 
Figs. 3,4). In its turn, the latter requirement is based on the theoretical results 
that r(fc, ±A) is a smooth monotonous function of k changing between and 1 (see 
Ref.^). Therefore in the present paper we restrict ourselves to considering elec- 
trons and positrons for M = 10^^ g. Including the heavier fundamental fermions 
(/U^-mesons and r^-leptons) within the given algorithm is only possible for smaller 
M. 

4. Estimations for luminosities 

Before presentation of the numerical results we should touch upon the conver- 
gence of the series (11)-(12) over I and n respectively. For this aim we denote 



Cl 



00 

(") = / 



|5ii(fc,A)p + |5ii(fc,-A)p 

e87rfe+ + 1 



■k+dk 



and represent the coefficients c;(n) in the form (omitting the integrand) 

(ln;)/27r 00 

ci{n) = cn{n) + ci2{n) = j ^ J ' 

M (lnZ)/27r 



(25) 
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SO that L{n) of (11) is equal to Li{n) + L2{n) respectively. Also it should be noted, 
as follows from the results of Ref.^, we have for barrier transparency r(A;,±A) = 
|sii(/c,±Ap that r(/i, ±A) = and at fc — > +oo 

r(/c, ±A) = |sn(fc, ±A)|2 = 1 + O(fc-i) . (26) 

Further we have the asymptotic behaviour for large I and k « I 



\s,,{k,±X)\=C (k-i_i)[l + o{l)] (27) 

^/n + / + 1 

with some constant C. To obtain it wc should work using the equations (18)-(21) 
by conventional methods of mathematical analysis so that the detailed derivation 
of (27) lies somewhat out of the scope of the present paper and we will not dwell 
upon it here. Using (25) and (27) we find 



-2Vn+i+l /•(ln/)/27i- _ ..\2i.+ 

c,,{n) < B- -— / ^\ J'J \ dk 

with some constant B and then it is clear that 



10 

+1 r° 

n + l + 1 Jo 6^'^'=+ - l"" (87r)4 n + l + 1 ' 

where we used the formula (for natural p > 0)^ 



^ ^^^2Vn+m TOO (fc+)3 ^^^^ _ ^3!C(4) e-2v'^^+^ 



OC 

/ 



p!C(p+l) 



e*- 1 





with the Riemann zeta function C(*)) while ({4) = Tr'^/90. As a consequence, one 
may consider 



~e-^(v5rTT--^), ,28) 
where we employed the asymptotical behaviour of the incomplete gamma function® 

CO 

r(a,a;) = ^ e'H'^-^dt ~ a^^-^e"^, a; ^ +oo , 

X 

SO that the series is evidently convergent. 
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At the same time due to (26) 

oo 

ci2{n)<D J e-^''''^k+dk = 

(lni)/27r 

^ 2 

D j e-8"'=^^^dfc+ , z = V(ln;/27r)2-/x2 

z 

with some constant D. Since < 1 and 



then 



This entails 



(lii0/27r 



cn{n) < D^^ . (29) 



C50 



, 2£+l)lnZ Inn 
^2(n) ~ / dl ~ ^ . (30) 



and the series X]i°-^2('^) is also convergent. Under this situation we obtain that 
the series of (12) is convergent, i. e., L < +oo. Consequently, we can say that all 
the further computed luminosities exist and are well defined. 

5. Numerical results 

In view of (26), we have 



T{k,X)k+dk r kdk 
e8^fe+ + l ~y gS^fe-l ~ (87r)2^'^' ~ 6-82 



C(2) = ^ , (31) 



whilst the latter integral can be accurately evaluated using the trapezium formula on 
the [0, 3] interval. Having confined the range of k to the [/x, 3] interval, accordingly, 
it is actually enough to restrict oneself to < /, n <15-20 when computing L{n), L. 
We took <l < 20, < n < 15. Since the Wronskian of (17) does not depend on x, 
the latter should be chosen in the region where the potentials of (22) are already 
small enough. Besides, potentials q~ of (22) are really equal to when x < -20. We 
computed r(fc. A) according to (17) B.t x = xq — 300, where were obtained from 
the Volterra integral equations, respectively, (18) and (23). For this aim, according 
to the methods developed for numerical solutions of integral equations (see e. g., 
Ref.^), those of (18) and (23) have been replaced by systems of linear algebraic 
equations which can be gained when calculating the conforming integrals by the 
trapezium formula, respectively, for the [-20, 300] and [300, 400] intervals. The 
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sought values of / were obtained as the solutions to the above linear systems. 
After this the derivatives (/^)^ were evaluated in accordance with (19) and (24) 
while employing the values obtained for The typical behaviour of r(fc, A) is 
presented in Figs. 3, 4. 

Figure 3: Typical behavior of barrier transparencies for A > 0. 



Figure 4: Typical behavior of barrier transparencies for A < 0. 

Further Fig. 5 presents luminosity Li,(0)/A with A of (11) for electronic neutrino 
(antineutrino) (no interaction with monopoles) as function of k and also the un- 
twisted Le{0)/A and the all configurations L^/A luminosities for electron (positron) 
as functions of k. The areas under the curves give the corresponding values of 
L{0)/A and L/A. 

Figure 5: Luminosity for electronic neutrino (antineutrino) and untwisted and all 
configurations luminosities for electron (positron). 

In Table 1 the data on computation of the untwisted and the all configurations 
luminosities are represented for all particle species under consideration. Finally, 
we can introduce the total untwisted Cq (as the sum over all species in the second 
column of Table 1) and the total C (as the sum over all species in the third column 
of Table 1) luminosities multiplied by A that should be detected by an external 
observer near black hole. We shall have 

Co = 0.100317 • 10^^ erg ■ s'^ C = 0.114389 • 10^*^ erg • s'^ 

so that contribution owing to Dirac monopoles amounts to 12.3022 % of £. 

6. Concluding remarks 

It seems to us the considerations of the present paper show that exact quan- 
titative study of the Hawking radiation is not so simple task. At the same time 
such an exploration is important enough in a number of astrophysical problems, 
for example, for more precise determination of the masses of primordial black holes 
which just expire today. But, as we have seen, just the black hole mass is one of 
the crucial parameters that define both the relevant scattering problems and the al- 
gorithms necessary for numerical calculations. Besides the exact evaluations might 
shed a light on nontrivial topological properties of black holes, in particular, on 
those connected with residence of Dirac monopoles (and generally U(N)-monopoles 
with > 1) in black holes. As a result, one should continue the exact calculations 
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Table 1. Untwisted and the all configurations luminosities. 



Particle 


L(i))/A 




L/A 




]\Ioiioi)ol(> (:oiitril)utioii of L (%) 




0.786745-10- 


-;i 


0.786745-10- 


-;5 







0.786745-10- 


-3 


0.786745-10- 


-3 





e~ 


0.120798-10- 


-2 


0.148780-10- 


-2 


18.8072 


e+ 


0.120798-10- 


-2 


0.148780-10- 


-2 


18.8072 



within the broader spectrum of both black hole masses and quantum particles emit- 
ting from black holes, in first turn, for other fundamental fermions - /it^-mesons and 
r^-leptons. Also one should pass on to the other types of black holes, for instance, 
to the Roissncr-Nordstrom ones where the description of the necessary S-matrices 
for spinor particles has been just recently obtained. 
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